Parameter Inference in the Pulmonary Circulation of Mice

This study focuses on parameter inference in a pulmonary blood cir- culation model for mice. It utilises a fluid dynamics network model that takes selected parameter values and aims to mimic features of the pulmonary haemody- namics under normal physiological and pathological conditions. This is of medical relevance as it allows monitoring of the progression of pulmonary hypertension. Constraint nonlinear optimization is successfully used to learn the parameter values

MCMC with Delayed Acceptance using a Surrogate Model with an Application to Cardiovascular Fluid Dynamics

Parameter estimation and uncertainty quantification in physiological modelling is a vital step towards personalised medicine. Current state-of-the-art in this research area performs parameter optimisation, with very limited uncertainty quantification. This paper demonstrates the advantage of novel sampling methods, applied on a complex biological PDE system of the pulmonary circulation. The aim is to find an efficient and accurate method for the inference and uncertainty quantification of unknown parameters, relevant for disease diagnosis (pulmonary hypertension) from limited and noisy blood pressure data. The data likelihood is expensive to evaluate as it requires solving numerically a system of PDEs. Therefore, having a model that best trades off accuracy and computational efficiency is of uppermost importance. In this study, we employ fast Bayesian methods, namely MCMC algorithms coupled with emulation using Gaussian Processes, to achieve a computational speed-up. We compare the Delayed Rejection Adaptive Metropolis algorithm in a History Matching framework to the delayed acceptance Adaptive Metropolis algorithm. The first algorithm draws samples from the approximate posterior distribution, while the latter is guaranteed to generate samples from the exact posterior distribution asymptotically. In this paper we propose and derive the n-steps ahead delayed acceptance Metropolis-Hastings algorithm, which is a generalisation of the classical 1-step ahead delayed acceptance Metropolis-Hastings. We show the superiority of the n-steps ahead algorithm over the 1-step ahead method

Numerical simulation of blood flow and pressure drop in the pulmonary arterial and venous circulation

A novel multiscale mathematical and computational model of the pulmonary circulation is presented and used to analyse both arterial and venous pressure and flow. This work is a major advance over previous studies by Olufsen et al. (Ann Biomed Eng 28:1281–1299, 2012) which only considered the arterial circulation. For the first three generations of vessels within the pulmonary circulation, geometry is specified from patient-specific measurements obtained using magnetic resonance imaging (MRI). Blood flow and pressure in the larger arteries and veins are predicted using a nonlinear, cross-sectional-area-averaged system of equations for a Newtonian fluid in an elastic tube. Inflow into the main pulmonary artery is obtained from MRI measurements, while pressure entering the left atrium from the main pulmonary vein is kept constant at the normal mean value of 2 mmHg. Each terminal vessel in the network of ‘large’ arteries is connected to its corresponding terminal vein via a network of vessels representing the vascular bed of smaller arteries and veins. We develop and implement an algorithm to calculate the admittance of each vascular bed, using bifurcating structured trees and recursion. The structured-tree models take into account the geometry and material properties of the ‘smaller’ arteries and veins of radii ≥ 50 μ m. We study the effects on flow and pressure associated with three classes of pulmonary hypertension expressed via stiffening of larger and smaller vessels, and vascular rarefaction. The results of simulating these pathological conditions are in agreement with clinical observations, showing that the model has potential for assisting with diagnosis and treatment for circulatory diseases within the lung

Influence of image segmentation on one-dimensional fluid dynamics predictions in the mouse pulmonary arteries

Computational fluid dynamics (CFD) models are emerging as tools for assisting in diagnostic assessment of cardiovascular disease. Recent advances in image segmentation has made subject-specific modelling of the cardiovascular system a feasible task, which is particularly important in the case of pulmonary hypertension (PH), which requires a combination of invasive and non-invasive procedures for diagnosis. Uncertainty in image segmentation can easily propagate to CFD model predictions, making uncertainty quantification crucial for subject-specific models. This study quantifies the variability of one-dimensional (1D) CFD predictions by propagating the uncertainty of network geometry and connectivity to blood pressure and flow predictions. We analyse multiple segmentations of an image of an excised mouse lung using different pre-segmentation parameters. A custom algorithm extracts vessel length, vessel radii, and network connectivity for each segmented pulmonary network. We quantify uncertainty in geometric features by constructing probability densities for vessel radius and length, and then sample from these distributions and propagate uncertainties of haemodynamic predictions using a 1D CFD model. Results show that variation in network connectivity is a larger contributor to haemodynamic uncertainty than vessel radius and length

A computational framework for generating patient-specific vascular models and assessing uncertainty from medical images

Patient-specific computational modeling is a popular, non-invasive method to answer medical questions. Medical images are used to extract geometric domains necessary to create these models, providing a predictive tool for clinicians. However, in vivo imaging is subject to uncertainty, impacting vessel dimensions essential to the mathematical modeling process. While there are numerous programs available to provide information about vessel length, radii, and position, there is currently no exact way to determine and calibrate these features. This raises the question, if we are building patient-specific models based on uncertain measurements, how accurate are the geometries we extract and how can we best represent a patient's vasculature? In this study, we develop a novel framework to determine vessel dimensions using change points. We explore the impact of uncertainty in the network extraction process on hemodynamics by varying vessel dimensions and segmenting the same images multiple times. Our analyses reveal that image segmentation, network size, and minor changes in radius and length have significant impacts on pressure and flow dynamics in rapidly branching structures and tapering vessels. Accordingly, we conclude that it is critical to understand how uncertainty in network geometry propagates to fluid dynamics, especially in clinical applications.Comment: 21 pages, 9 figure

Representation of Time-Varying Stimuli by a Network Exhibiting Oscillations on a Faster Time Scale

Sensory processing is associated with gamma frequency oscillations (30–80 Hz) in sensory cortices. This raises the question whether gamma oscillations can be directly involved in the representation of time-varying stimuli, including stimuli whose time scale is longer than a gamma cycle. We are interested in the ability of the system to reliably distinguish different stimuli while being robust to stimulus variations such as uniform time-warp. We address this issue with a dynamical model of spiking neurons and study the response to an asymmetric sawtooth input current over a range of shape parameters. These parameters describe how fast the input current rises and falls in time. Our network consists of inhibitory and excitatory populations that are sufficient for generating oscillations in the gamma range. The oscillations period is about one-third of the stimulus duration. Embedded in this network is a subpopulation of excitatory cells that respond to the sawtooth stimulus and a subpopulation of cells that respond to an onset cue. The intrinsic gamma oscillations generate a temporally sparse code for the external stimuli. In this code, an excitatory cell may fire a single spike during a gamma cycle, depending on its tuning properties and on the temporal structure of the specific input; the identity of the stimulus is coded by the list of excitatory cells that fire during each cycle. We quantify the properties of this representation in a series of simulations and show that the sparseness of the code makes it robust to uniform warping of the time scale. We find that resetting of the oscillation phase at stimulus onset is important for a reliable representation of the stimulus and that there is a tradeoff between the resolution of the neural representation of the stimulus and robustness to time-warp. Author Summary Sensory processing of time-varying stimuli, such as speech, is associated with high-frequency oscillatory cortical activity, the functional significance of which is still unknown. One possibility is that the oscillations are part of a stimulus-encoding mechanism. Here, we investigate a computational model of such a mechanism, a spiking neuronal network whose intrinsic oscillations interact with external input (waveforms simulating short speech segments in a single acoustic frequency band) to encode stimuli that extend over a time interval longer than the oscillation's period. The network implements a temporally sparse encoding, whose robustness to time warping and neuronal noise we quantify. To our knowledge, this study is the first to demonstrate that a biophysically plausible model of oscillations occurring in the processing of auditory input may generate a representation of signals that span multiple oscillation cycles.National Science Foundation (DMS-0211505); Burroughs Wellcome Fund; U.S. Air Force Office of Scientific Researc

On the Mechanics Underlying the Reservoir-Excess Separation in Systemic Arteries and their Implications for Pulse Wave Analysis

Several works have separated the pressure waveform p in systemic arteries into reservoir pr and excess pexc components, p = pr + pexc, to improve pulse wave analysis, using windkessel models to calculate the reservoir pressure. However, the mechanics underlying this separation and the physical meaning of pr and pexc have not yet been established. They are studied here using the time-domain, inviscid and linear one-dimensional (1-D) equations of blood flow in elastic vessels. Solution of these equations in a distributed model of the 55 larger human arteries shows that pr calculated using a two-element windkessel model is space-independent and well approximated by the compliance-weighted space-average pressure of the arterial network. When arterial junctions are well-matched for the propagation of forward-travelling waves, pr calculated using a three-element windkessel model is space-dependent in systole and early diastole and is made of all the reflected waves originated at the terminal (peripheral) reflection sites, whereas pexc is the sum of the rest of the waves, which are obtained by propagating the left ventricular flow ejection without any peripheral reflection. In addition, new definitions of the reservoir and excess pressures from simultaneous pressure and flow measurements at an arbitrary location are proposed here. They provide valuable information for pulse wave analysis and overcome the limitations of the current two- and three-element windkessel models to calculate pr